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THE APPROXIMATION OF STELLAR ENERGY 
DISTRIBUTIONS AND MAGNITUDES FROM 
MULTI-COLOR PHOTOMETRY 

By 

Paul B. Davenport 


SUMMARY 

A procedure is developed for approximating the energy distribution 
of a star from present day multi-color photometry. This approximation 
is then used to predict stellar magnitude relative to a light sensitive 
device with known spectral response. Using magnitudes in the U, B, V 
system it is found that the energy distributions of stars can be approxi- 
mated by this procedure fairly accurately over the range of 3000-8000 
angstroms. The stellar magnitude relative to a sensor whose dominant 
response lies in this region can then be determined with an accuracy of 
about 0.1 magnitude. 
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THE APPROXIMATION OF STELLAR ENERGY 
DISTRIBUTIONS AND MAGNITUDES FROM 
MULTI-COLOR PHOTOMETRY 


INTRODUCTION 

With the increasing applications of star trackers in air and space navigation 
it is imperative that procedures be developed which accurately predict the magni- 
tude of a star relative to a given tracker. It certainly is not practical, nor in 
general feasible, to observe every star of interest each time a tracker with a 
different response becomes available. There have been many internally published 
documents from various organizations giving the photoelectric magnitudes of the 
brighter stars. All of these lists, however, have been generated in practically the 
same manner i.e. black body radiation was assumed (the temperature obtained 
from various sources as a function of the Draper spectral type) to generate an 
index which is added to the Harvard visual magnitudes. Unfortuntely, this assump- 
tion as well as the basic data used is not accurate when compared with modern 
observations. These methods have also assumed that all sensors with the same 
S-number have the same spectral response which can be obtained from a hand- 
book containing the nominal response. The fact is, however, that even when two 
devices from the same manufactor are designated the same their measured re- 
sponse may be quite different causing deviations of up to a tenth of a magnitude. 
(This value was derived from our own calculations based upon the measured 
responses of dozens of photomultipliers of the same type. See also [1]). 

Here, we present a procedure for determining stellar magnitudes which uses 
the latest available data, namely multi-color photometry such as the U, B, V 
system [2], [3 ], and [4 ]. The approach here is to obtain information about the 
energy distribution of a star from the observed color magnitudes which in turn 
can be used to determine the stellar magnitude relative to a sensor whose re- 
sponse lies in the region covered by the responses defining the color magnitudes. 
Since a large number of stars have now been observed in the U, B, V system, 
and the dominant part of the response curves of many sensors lies in that part 
of the spectrum covered by the U, B, and V filters the techniques developed here 
have a very general use. 


THE MAGNITUDE EQUATION 

The stellar magnitude (m s ) of an astronomical radiation source reduced to 
outside the earth's atmosphere and incorporating the interstellar absorption is 
given by 
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m s0 - 2.5 log—— = “ 2.5 log I s + k s , 

J- rt 


J(M o-(K) d\ . 


J(A) is the spectral distribution of light intensity from the source, a s is the 
spectral response function of the optical train observing the radiation, and the 
constant I . is the value of I for a standard source of radiation with magnitude 

S 0 s 

m sQ . If a black -body distribution is assumed then 


jw = 


: 2/T\ _i 


where T is the temperature of the source. 

To obtain information about the spectral distribution of light intensity from 
a multi-color magnitude system such as Johnson's & Morgan's UBV system we 
use the constant energy concept widely practiced in photometric work. Since 
many observers today use the effective wave number concept we adopt this pro- 
cedure here [5]. 

Let x = 1A , then the effective intensity integral equation (2) becomes 


X -2 J (x) CT (x) dx 


‘'a 


) <lX x ) dx 


where 


F(x) = x“ 2 J(x) . 
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For a black-body distribution 


F(x) = _fi X 


, c 2 /T X _i 


(4) 


Letting 


x o = 


J roo 

xcr 

0 * 


(x) dx 


f 

Jo 


cr (x) dx 


and expanding F(x) in a Taylor series about x Q the effective intensity integral 
becomes 


I = 


( x o> r 

^0 


F ( x o> I cr s ( x ) dx+ ^~ 
2 dx 2 


f 


( x ~ x 0 ) 2 °" s ( x ) dx + 


x = x 0 ^0 


f o- s (x) dx [f(x 0 > + Af] 
Jo 


where 


AF 


L I d n F 
dx n 


n “ 2 


^n( X 0> 


and 


M n ( x 0 ) = 


J roo 

1 ( x ■ x 0 ) n 
0 


cr ( X ) dx 


J 0 a . (x) dx 


The stellar magnitude is then given by 


m 


= -2.5 log [f(x 0 ) +AF] +k' s , 


= -2.5 [log F(x 0 ) + log (1 + e)] +k' s , 
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where 


( 5 ) 


_ AF 
£= F(X 0 ) 

Since 

log (1 + e) = (log e) ~\ &2 + • • > 

provided | e | < 1 


m = - 2.5 


log F(x 0 ) + (log e) 



+ k' 


Thus, if the band pass of the receiver is narrow enough or the second and 
higher derivatives of the energy distribution F(x) are sufficiently small the 
magnitude is a measure of the monochromatic flux at x Q i.e. 


m s = - 2.5 log F(x 0 ) + k' s • 

For the U, B, and V bandpasses we assume that this is the case. Hence, 


U = - 2.5 log F(Xj) + k x 

B = - 2.5 log F(x 2 ) + k 2 (6) 


where 


V = - 2.5 log F(x 3 ) + k 3 


X. 

1 



Xcx.(X) dx 



i = 1. 2, 3 


(7) 


and o\ (X) is the corresponding response function. Thus the U, B, V system 
provides a means of approximating the energy distribution at three points within 
the response of many photomultipliers and other sensors. 

Since the band width of most sensors is at least three times as wide as each 
filter (U, B, or V) we assume that the energy distribution over the range of 
interest is quadratic. 
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Thus, the S magnitude, m s , is given by 

m s = - 2. 5 [log F(X s ) + (log e) e] + k s (8) 

where e is defined by equation (5) with the third and higher derivatives assumed 
to be zero, x s is given by equation (7) and the response function used in the 
equations is that of the measured response of the sensor of interest. 

From the three values of y(x) = log F(x) obtained by the U, B, and V magni- 
tudes (equation 6) a quadratic approximation to y(x) can be determined by con- 
structing a quadratic function passing through the three points i.e. 

y(X) = AjX 2 + a 2 x + a 3 , (9) 

where the coefficients A. are determined from the matrix equation 

/*?■*! A /a A y( x i)\ fy\ 

X 2 X, 1 A, = y(x 2 ) = y 2 

\ a 3/ \y( x 3>/ \y 3 / . 

The explicit solution is given by 

A i = - [ (x - x 3 Xy 3 - y 2 ) + ( x 3 “ x 2 )(y 3 - y»)] - 

A 2 = [( x 3 ~ x i )(y 3 ■ y 2 ) + ( x 2 ‘ x 3>(y3 ~ yi) > ( 10 ) 

A 3 = ^ [ x 2 x 3 ( x 2 " x 3>yi + x i x 3 ( x 3 - x i)y 2 + x i x 2 ( x i - x 2 )y 3 ] - 

where 

D = (X 3 - X 2 )(X 3 - X 1 )(X 1 - X 2 ) . 
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To determine the S magnitude as given by equation (8) we need y(X s ) = logF(X s ) 
and e. y(X s ) is obtained directly from equation (9) i.e. 


Now 


where 


Since 


y(X s ) = log F(X s ) = A x X 2 + a 2 x s + a 3 
1 d 2 F 


2 dX 2 


e =• 


AF 

F(XJ 


F(X S ) 


M 2 (X s ) 


^ 2 ( X s> = 


r 

Jo 


(X - X ) 2 CT (X) dx 


/*00 

a(X) dx 

Jo 


F(X) = 10 y ( x) , 


(11) 


then 


— = -In 10 10 y < x > ^ , 
dX dX 


^ 10 
dX 2 


d 2 F 
*2 


— 

/ \2 

l 0 y(x) 10 10 y ( x > 

( — ] 

dX 2 

\ dX / 


- = 'in 10 
F (X) 


+^niofeV 
dx 2 W 


From equation (9) 


2 A, X + A, , 

Is 2 . 


thus 


6 = — In 10 
2 


dy 
dX 

dV 

dX 2 

2A X + In 10 (2AjX s + A 2 ) 2 


= 2A, 


/- 2 ( X s) 


( 12 ) 
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In summary, the constants X., X,, X,, and X are determined from the 
known response of the U, B, V, and S bandpasses respectively (^ 2 (X 2 ) is 
obtained from the S response and x s ) . y x , y 2 , and y 3 are obtained from equation 
(6), then A r A 2 , and A 3 from equation (10) log F(X s )from equation (11), € from 
equation (12), and finally the S magnitude (m s ) from equation (8). 

By convention the values of U, B, and V are not given directly, but as a single 
magnitude (V) and two color indices (B-V and U-V or U-B). Since the data is 
often obtained in terms of V, B-V, and U-V [4] we will write our solution in terms 
of these variables. Thus 

A, = 215 [ <X > " XjKB ' v - V * < X J * V < u ' v - c „> 


a 2 = Tsv [< x ! - x i)< B • v - C B> + (Xf - x i>( u - v - c„>] 

A, = 2^ j (k - ' V) ~d ■ X 3 )<U • V " C -> 

+ X 1 X 3 (X 3 -X 1 )(B-V-C b )]| 

where 

C u = k i - k 3 , Cg = k 2 - k 3 • 


Therefore, 

m s = V + b x (U - V) + b 2 (B - V) + b 0 + Am , (13) 


where 

xl - X S (X 2 + X 3 ) + X 2 X 3 + fi 2 (X s ) 
1 = (X x - X 2 )(X 1 - x 3 ) 


Xg - X,(X 1 + x 3 ) +x t x 3 + /X 2 (X ) 
(X 3 - XaXX, - XJ 


Am = - -i 'In 10 yU- 2 (X s ) 


C^U - V) + C 2 (B ~ V) + c 3 
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X, - 2X + X, 

1 S <3 


C. = 


X„ - 2X + X. 


s 3 r _ 

. ^2 “ 


1 -(X 1 -X 2 )(X 1 -X 3 ) ’ 2 (X 1 -X 2 )(X 3 -X 2 ) 


C, = (CCl + CgC 2 ) 


The constant b 0 is arbitrary and serves as a zero point for the magnitude 
system. The only constant dependent upon stellar energy distributions is C 3 
which does not require absolute measures but only relative ones i.e. ultraviolet 
and blue relative to visual. Besides, the constant C 3 is involved only in the second 
order term. Since the effective wave numbers are also relative measurements 
it appears that the constants in equation (13) can be determined quite accurately. 
The primary source of error, therefore, will probably be due to the assumptions 
that the flux is linear over each bandpass and quadratic over the entire interval 
of interest. 


CONSTANT ENERGY EFFECTIVE WAVE NUMBERS 

Using the published response curves of the U, B, and V filters [ 2 ] and the 
reflectivity of two aluminum reflections [ 6 ] we obtain 

Xj = 2.889, Xj = 2.296, x 3 = 1.826, 

which is consistent with Code's results L5 ]. 

If the energy distributions were linear over each bandpass then the plots of 
the distributions versus wave number for different stars (with the same magni- 
tude) would intersect at a point. This point of intersection would give the effective 
wave number and the value of the flux at the effective wave number. Actually 
the plots will not intersect at a point since the flux is not linear, however, a point 
of minimum deviation will exist which approximates the effective wave number. 

Using Willstrop's data for V=0.0 [73 we find that the intersection occurs near 
X= 1.83 andF(1.83) = 1.14X 10~ 5 (erg/cm 2 /sec) /j. . The effective wave number ob- 
tained on this manner agrees with that obtained directly from the V sensitivity- 
curve. The value of F(1.83) is consistent with Code's value of 1.13 x 10 5 (erg/ 
cm 2 /sec) m using Minnaert's derived fluxes and Stebbins & Kron's apparent visual 
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magnitude of the sun. It is also consistent with Code's direct measurement of 
1.20 at 1.80. We shall adopt the value of 1.14 x 10 -5 (erg/cm 2 /sec) /x at 1.83 
(l//u) for V= 0.0. Thus, 

V = - 2.5 log F(1.83) + 0.14 


From Willstrop's data for B = 0.0 the effective wave number of the B band- 
pass is 2.295 and the average value of F(2.295) is 1.24 x 10" 5 (erg/cm 2 /sec)M. 
Using the adopted expression for V and Code's expression for the monochromatic 
magnitudes we obtain 


log F(X) = 0.4 |m(1.83) - m(X) - V + 0.142J 

where m(x) is Code's monochromatic magnitudes. From Code's data using the 
observed values of B-V and selecting V such that B = 0.0 we find that the effective 
wave number for the B band pass is 2.295 and the average value of F(2.295) is 
1.26 x 10" 5 (erg/cm 2 /sec) M for B = 0.0. Oke's data [8 .] is for the stars whose 
spectral types are very similiar hence the plotting method using this data does 
not yield a satisfactory approximation to the effective wave numbers i.e. the 
curves are nearly parallel. The average value of the flux at the effective wave 
numbers, however, should be consistent with that of Willstrop's and Code's. 

Indeed this is the case, using Oke's monochromatic magnitudes in the same manner 
as Code's we obtain an average value of 1.29 x 10 “ 5 (erg/cm 2 /sec ) /x at 2.29 (1/m) 
for B = 0.0. Thus, the effective wave number as obtained from Johnson's, Willstrop's, 
and Code's data is in agreement and the value for the flux at the effective wave 
number from Willstrop's, Code's, and Oke's data is also consistent. Hence, 

B = - 2.5 log F( 2. 295) + 0.251 . 

From Code's data, the observed values of U-V, and selecting V such that 
V = 0.0 we obtain an effective wave number for the U band pass of 2.905 with an 
average value of 0.50 x 10” 5 (erg/cm 2 /sec) M for the flux. Oke's data gives an 
average value of 0.41 (erg/cm 2 /sec) m at 2.90. Thus, we adopted the value of 
0.45 (erg/cm 2 /sec) /x at 2.89 (1 /a 0 for U = 0.0. Hence, 


2.5 log F(2.89) - 0.867 


A summary of the results of the various investigators and the adopted values 
are given below. 
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Effective Wave Number 


u 

B 

V 

Remarks 

2.89 

2.29 

1.83 

Sensitivity — curves by Johnson & 
Morgan — calculated by Code 

2.889 

2.296 

1.826 

Sensitivity — curves by Johnson & 
Morgan - calculated by author 


2.295 

1.83 

Willstrop's data 

2.905 

2.295 


Code's data 

2.89 

2.295 

1.83 

Adopted 


Flux at Effective Wave Number 

U = 0 
X = 2.90 

B = 0 
X = 2.295 

V = 0 
X = 1.83 




1.13 

Code using Minnaert's flux and 
Stebbins & Kron V of sun 


1.24 

1.14 

Willstrop's data 

0.50 

1.26 

1 


Code's data 

0.41 

1.29 


Oke's data 

0.45 

1.26 

1.14 

Adopted 


AN EXAMPLE 

From the nominal S4 response [ 9 ] we obtain as the effective number of the 
S4 response the value 2.49 and m- 2 (2.49) = 0.155. With these values plus the 
values given above we obtain 

b x = 0.439, b 2 = 0.413, 

Cj = - 1.321, c 2 = 0.890, 

c 3 = - 1.43 
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Hence, 


m s4 = V + 0.439(11 - V) + 0.413(B - V) + 0.145 (14) 


-0.071 [- 1.32(U - V) + 0.89(B - V) - 1.43] 2 

where the zero point has been selected such that m s4 - V = 0 when U-V = B- V=0 

To test the accuracy of this formula we generated the flux of seven stars 
(e Ori, /3 Ori, a Lyr, /3 Ari, a Boo, X Ser & a Tau) from Code’s data. The seven 
stars were selected so as to cover a wide range of U - V and B -V values. The values 
of V used to generate these data were chosen so that B = 0 merely to keep the 
values of the flux within the same numerical range. 

The values beyond X = 2.94 were obtained by linear extrapolation using the 
values at X = 2.74 and X = 2.94. Values at even increments were obtained by 
linear interpolation. Using these energy distributions and the response curves the 
deflection y, the blue-yellow color index c b , and the ultraviolet-yellow color 
index were computed using numerical integration. Values of V, B - V, and 
U - V were then computed using a zero point so as to minimize the mean differ- 
ence between the observed and computed values. 

The relationships were 


V = - 2.5 logy + 0.22 , 


U - V = 2.5 log-^-- 0.05 , 

B - V = 2.5 log^- + 0.79 . 
b 


Except for the U - V of /3 Ori and a Boo the agreement between the observed 
and computed values are quite good and the extrapolation in the ultraviolet can 
easily account for the error of 0.1 mag. in U - V for these two stars. 

The energy distributions from Code were then assumed to be correct and the 
computed V, U - V, and B - V (from these energy distributions plus the U, B, and 
V response functions) to be correct also. The S4 magnitudes were then computed 
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by integration and by equation (14), (actually the value of C and the zero point 
were adjusted to agree with Code’s data rather than the adopted means). The 
maximum error is only 0.06 magnitudes which includes uncertainties in the 
sensitivity-curves as well as the assumed energy distributions. Since the con- 
stants used in equation (14) are based on data which is consistent from several 
different sources it appears that the S4 magnitude can be predicted from the 
observed V, U - V, B - V magnitude & colors to within 0.1 magnitude (assuming 
that the S4 spectral response is accurately known). 

Actually, the energy distributions can be approximated much better by using 
the relationships 


F(X) = 10 y < x > , 


y(x) = AjX 2 + A 2 X + A 3 , 


rather than retaining only quadratic terms in the expansion of 1 0 y ' x ^ . This 
additional accuracy in the energy distrituions, however, affects the S4 magnitudes 
very little. This is due to the fact that the third and higher moments (/x 3 ,/x 4 . . . ,) 
about the effective wave number for this particular curve are very small. For a 
response curve with a wider band pass of a large third moment about the effective 
wave number the better approximation may be necessary. 


CONCLUSION 

The method presented here was specifically formulated with the intent of 
being applied to a S4 response and particularly to obtain the guide star list for 
the Orbiting Astronomical Observatory launched in April of 1966. As our example 
shows, the method is quite accurate for this type of response. Since that time 
we have found that as suspected the accuracy falls off for sensors which respond 
to a wider range of light. However, with more recent data and minor modifications 
to our method here it is reasonable to expect that similar procedures will be 
applicable for any sensor whose response lies in the 3000-8000 Angstrom range 
with an accuracy of better than 0.1 magnitudes. The lesults of these late improve- 
ments will be presented at a later date. 
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